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ABSTRACT 

Observations with the Swift satellite of X-ray afterglows of more than a hun- 
dred gamma ray bursts (GRBs) with known redshift reveal ubiquitous soft X-ray 
absorption. The directly measured optical depth r at a given observed energy 
is found to be constant on average at redshift z > 2, i.e., (r(0.5 keV)) 2>2 = 
0.40 ± 0.02. Such an asymptotic optical depth is expected if the foreground 
diffuse intergalactic medium (IGM) dominates the absorption effect, and if the 
metallicity of the diffuse IGM reaches ~ 0.2 — 0.4 solar at z — 0. To further test 
the IGM absorption hypothesis, we analyze the 12 highest S/N (> 5000 photon) 
z > 2 quasar spectra from the XMM-Newton archive, which are all extremely ra- 
dio loud (RLQs). The quasar optical depths are found to be consistent with the 
mean GRB value. The four lowest-z quasars (2 < z < 2.5), however, do not show 
significant absorption. The best X-ray spectra of radio-quiet quasars (RQQs) at 
z > 2 provide only upper limits to the absorption, which are still consistent with 
the RLQs, albeit with much lower S/N (<1000 photons at z w 4). Lack of quasar 
absorption poses a challenge to the smooth IGM interpretation, and could allude 
to the opacity being rather due to the jets in RLQs and GRBs. However, the jet 
absorbing column would need to appear in RLQs only at z > 2.5, and in GRBs to 
strongly increase with z in order to produce the observed tendency to a constant 
mean r. High X-ray spectral resolution can differentiate between an absorber 
intrinsic to the source that produces discernible spectral lines, and the diffuse 
IGM that produces significant absorption, but no discrete features. 

Subject headings: Cosmology, 98.80.-k, Gamma ray bursts, 98.70.Rz, Quasars, 98.54.Aj 



1. Introduction 

Swift X-ray spectra of gamma ray bursts (GRB) afterglows reveal prevalent soft X- 
ray absorption, which is commonly assumed to originate in the host galaxy. The standard 
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X-ray absorption measurement technique probes metal absorption, but quotes equivalent 
hydrogen column densities Nh by assuming a neutral absorber at the host (z), and with solar 
abundances. Under these assumptions, Njj towards GRBs shows a strong correlation with 
the host redshift z. The typical absorbing colu mn rises from Ny^ ~ 10 21 cm" 2 at z < 1 up to 



« 10 23 cm 2 at the highest observed redshifts ( 


Jako 


Dsson 


Watson et al. 


2007; 


Campana et al. 


2010; 


Rau et al. 


2010) 



20061 : ICampana et all 12006 



20101 ). Damped Lyman-a absorption is 



also seen in some GRB afterglows. Unlike the X-ray absorber, the redshift of the Lyman-a 
absorber is well c onstrained. The imp lied column is usually still well below the X-ray derived 
column (see, e.g.. Watson et al.ll2007l ). If the abundances in the X-ray absorber are sub-solar, 
the X-ray derived Ng values are even higher and the discrepancy with the Lyman-a column 
grows accordingly. Although an appreci able ionization rang e in a single medium can possibly 
account for such a discrepancy (recently. ISchady et al.ll2010l ). there is also the possibility that 
the X-ray and Lyman-a absorbers are physically distinct. 



High-z quasars al so commonly reveal soft X-ray absorption (e.g.. iFabian et al.l 12 001 



Worslev et al.ll2004al ]bl; IPage et al.ll2005l : lYuan et al.l200a iGrupe et al. 



2006 



Sambruna et al. 



20071 ). The implied absoring column, if intrinsic, is also of the order of 10 23 cm -2 . In the case 
of quasars as well, metal and hydrogen line absorption in the UV imply significantly lower 
Nh columns. In the case of low-luminosity active galaxies, partially ioni zed outflows are 
know n to have more X-ray column with only a trace of UV absorbing ions (ICrenshaw et al. 
20031 ). a discrepancy that again can be partially reconciled with an ionization corr ection. The 
X-ray absorption profile in three different quasars at z ~ 4.3 — 4.7 was noted by lYuan et al. 
(120051 ) to be remarkably similar, while ionized quasar outflows are not necessarily expected 
to be so uniform. Indeed, partially ionized outflows have not been directly identified in such 
luminous quasars, and thus it is plausible that the X-ray and UV absorbers towards high-z 
quasars are also physically distinct. 

Intrigued by the aforementioned puzzles associated with soft X-ray absorption of high- z 
sources, we wish to explore the similarities of X-ray absorption of high-z GRBs and quasars 
and their possible origi n. In Sec. [2] we present the X-ray opaci ties towards GRB afterglows 



from the Swift sample (jEvans et al. 



20091 ; ICampana et al.l 120101 ) . but now without assuming 



neither a redshift, nor an ionization, nor a metallicity for the X-ray absorber. We then discuss 
in Sec. [3] the viability of X-ray absorption by the diffuse intergalactic medium (IGM). In 
Sec. HJ we present a comparison sample of the highest signal to noise ratio (S/N) z > 2 quasar 
spectra drawn from the XMM-Newton archive. Realizing that our selection criterion based on 
the number of detected photons allows only radio loud quasars into the sa mple, in Sec. 14.1 
we also present th e best S/N high-z radio quiet quasars from the works of IShemmer et al. 
(120081 . |2006|, |2005j) . In Sec. [5] we conclude and compare the expected absorption signature in 
the X-ray spectra of high-z sources of a diffuse IGM versus that of jets and propose future 
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tests that can reveal whether the observed absorption is indeed due to the IGM or intrinsic 
to the sources. 



2. Observed-frame opacity of GRBs 



Soft X-ray absorption below 1 keV is mostly due to photo-ionization of heavy elements 
(e.g., C, N, O, Fe) and has been measured with the Swif t X-ray Telescope (XRT) in over a 



hund red GRBs with an identified redshift ( IGehrelsl 12004 ; Evans et al.ll2009t ICampana et al. 



20101 ) . It is preferable to use late photon arrival times when spectral variability is minimal 
in order to obtain the most reliable absorption measurements. In this work, we therefore use 
only the photon co unting (PC) mode data available from the Swift/GRB spectrum repository 
( Eyans et al.1 2009 , ]]) . T his choice of data is not as sophisticated as the careful time cuts 
of ICampana et al.l ( 120101 ). but on the other hand, it provides a uniformly reduced sample. 
We include all of the bursts up to 2010 July. We use only the 144 that have an identified 
redshift, out of a total of 520 GRBs. The sample greatly varies in S/N, but in order not to 
introduce any biases, we do not exclude any objects from the analysis. Out of the sample of 
144 GRBs, absorption at the 90% confidence level is measured for 113, while the other 31 
have only upper limits. 

The column density Njj(z) in the repository is conveniently given as an equivalent 
hydrogen column density assuming a neutral solar-abundance absorber at t he host redshift. 



The s trong correlation of these absorbing columns with z was noted by ICampana et al. 
(120101 ) . Here, we wish to carefully study the observed X-ray absorption effect, relaxing these 



assumptions as much as possible. We thus convert the quoted column into an optical depth 
at the observed energy E using 



r*(E) = a[(l + z)E]N H (z) 



where <r([l + z]E) is the total photo-ionization cross section per hydrogen atom at a photon 
energy of [1 + z\E and for solar metallicity gas. In Eq. [TJ we use the exact same cross section 
used by the reposit ory team, namely that prov ided by the phabs model in Xspe@ for a neutral 
solar-composition ( lAnders and Grevessd 119891 ) absorber. Importantly, the intrinsic column 
densities provided in the repository are th ose in excess of the no minal local Galactic column 
N^ al based on HI 21 cm measurements (jKalberia et al.l 120051 ). which is also given there. 



1 http:/ /www. swift. ac.uk/xrt_spectra/ 
2 http://heasarc. gsfc.nasa.gov/docs/xanadu/xspec/ 
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However, we note that the use of the solar abundances of lAnders and Grevessd (119891 ) for 
the Galactic column leads to an overestimate of the local X-ray absorption effect and, thus, to 
an underestimate of the residual (GRB) absorption, as these abunda nces are somewhat high 



compared to more up to date estimates of interstellar abundances dWilms et al.l l200oh . At 



0.5 ke V, e.g., the photo-ionization cross section per H atom based on lAnders and Grevesse 
( 119891 ) abundances wo uld be 7.14xl ~ 22 cm 2 compared to the preferred value of 6.22 xl0~ 22 
cm 2 obtained with the I Wilms et al.l (120001 ) abundances. This difference is demonstrated in 
Fig. [TJ We therefore introduce a small correction (increase) to the optical depth in Eq. [T] as 
follows: 



t{E) = r*(E) + N^ l [a A ^(E) - a w (E)] (2) 
where a A G (E) —a^{E) is the (positive ) differenc e in the inferr e d cros s section deduced from 



using the lAnders and Grevesse! ([1989, AG) and I Wilms et al.l (|2000| . W) abundances. On 



average, this correction to r is ~ 10%, but can reach 50% for high Galactic column GRBs. 

In Fig.[2]we present the PC-mean results for the optical depth r(0.5 keV) at the observed 
energy of 0.5 keV, where photo-electric absorption is significant and the instrument response 
is sufficiently high. The errors on r correspond to the 90% confidence errors quoted in the 
repository. For bursts whose absorption measurements are consistent with zero to within 
these errors, we plot only the +90% confidence (upper) limit. We stress that the optical 
depth is the directly measurable model-independent quantity, unlike column density. The 
scatter in r(0.5keV) in Fig. [2] is substantial, but there is no evolution with z, especially not 
at high z. Indeed, the error-weighted mean optical depth in Az = 1 bins is also shown in the 
figure, and its constancy with z is strikingly tight. For z > 2 bursts, the error-weighted mean 
optical depth is (r(0.5 keV)) 2>2 = 0.40 ± 0.02. The upper limits can not be included in an 
error-weighted mean, but a simple mean with no weighing yields (r(0.5 keV)) 2> 2 = 0.7±0.1 
with the upper limits, and (r(0.5 keV)) z> 2 = 0.55 ± 0.09 without them; Quoted errors are 
standard errors of the mean. Using the late-time PC spectra from the repository instead 
of the mean PC data increases the unweighed mean by ~ 0.1. This could be a result of 
the 2-3% steeper spectral slopes in the late-time spectra. On the other hand, excluding 
GRB 060202, a rare outlier with r(0.5keV) > 5, which also has a very steep spectral slope 
of r = 2.7 and is not fitted by the model very well, reduces the unweighed mean by ~ 0.1. 
All of these give a general idea of the systematic uncertainties that can be expected from 
the X-ray absorption measurements of the transient GRB afterglows. Given the wide range 
in the quality of the measurements, their detection limits, and the fact that including upper 
limits tends to increase the mean (i.e., largely inadequate measurements rather than low 
absorption), we prefer the error weighted mean of (r(0.5 keV)) 2 >2 = 0.40 ± 0.02, keeping in 
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mind, however, the large (~ 0.2) systematic errors. 

The unambiguous opacity at high- z and its apparent independence of z raises the sus- 
picion that the absorption is not intrinsic to the GRB, but is due to a rather uniform 
foreground. We focus on the high-;? population where host absorption is strongly quenched 
(as a sharply declines with E, see Sec. |3]). Therefore, at such high z the possibility that 
the GRB environment or host galaxy contribute to the observed absorption is small, unless 
the intrinsic column is of the order of Njj ~ io 22-23 cm -2 , well above the typical galactic 
column in the local universe. In order to illustrate the similarities of the absorption spectra 
irrespective of z, in Fig. [3] we plot spectral ratio plots for the twelve highest-z (3.85 - 8.26) 
GRBs with absorption confirmed at the 90% level. This sub-sample is drawn from the full 
Swift/XKT sample shown in Fig. [21 but is not necessarily representative of it. For Fig. El 



each spectrum is fitted to an absorbed power law (with IWilms et al.l I200CH . abundances), 
after which the excess (extra-galactic) absorption is removed. The ratio of the model (with 
only Galactic absorption) to the data, i.e., transmission, is plotted. The apparent drop of 
the ratio at low energies (E < 1 keV) thus reflects the additional photoelectric non-Galactic 
absorption towards each GRB. GRB 060202 is exceptional in Fig. [3] with its excessive ab- 
sorption, and can be seen also in Fig. [2] at 1 + z = 5.05 far above the mean GRB opacity. 
On the other hand, GRB 081029 at 1 + z = 4.85 is less absorbed (r « 0.16l{j;£|), although, 
formally, not as well constrained. For the most part, absorption sets in at E < 1 keV, and 
the transmission at E — 0.5 keV reaches ~ 0.6 — 0.7, i.e. r ~ 0.3 — 0.6. The column densities 
required to produce an opacity r that is fixed with z increase roughly as (1 + z) 2,5 due to 
the corresponding decrease of the photo-ionization cross section with energy, as can be seen 
in Fig. [1] and as explained i n mo re detail below. These columns then exceed 10 23 cm -2 as re- 



ported by lCampana et al.l ( 120101 ). In the next section, we explore an alternative explanation 



for the observed constancy of optical depth with z. 



3. Soft X-ray opacity of the diffuse IGM 

A natural origin of universal, isotropic, X-ray opacity that saturates at high-z is the 
diffuse intergalactic medium (IGM). In this section, we describe a simple diffuse IGM model, 
based on well established cosmological parameters, that can explain a tendency to a constant 
X-ray opacity for high- z sources. Since the photo-ionization cross section per H atom in the 
photon energy range of 0.5 keV < E < 10 keV scales roughly as cr{E) oc E~ 2 ' 5 (Fig. [TJ), 
for a redshifted absorber at a fixed observed energy E,a(E,z) oc (1 + z) -2 ' 5 . Also, since 
the absorption at X-ray energies is dominated by metals, one can assume cr(E, z, Z Q ) oc 
Z Q a(E, z) scales approximately with the IGM metallicity Z Q (in solar units) that can evolve 



-6 - 



with redshift as Z Q (z) = ZqT](z). We can thus write for the IGM photo-electric optical depth 



r IGM (E, z, Z Q ) = J n H {z')cr(E, z', Z Q )c \^pj dz' 



n ° cZ °a(E, 0) f (3) 



H ' J Jo (1 + z>)^{\ + z')sf{l + z'f Q M + n A 

where uq and Zq are the z = mean IGM hydrogen number density and metallicity (in solar 
units), c is the speed of light, and H is the Hubble constant. Q M and fl A are, respectively, the 
present-day matter and dark energy fractions of the critical energy density of the Universe. 
Even with no metallicity evolution, i.e., 7](z) = 1, the differential optical depth dr/dz' — > 
(1 + z')~ 2 at z 3> 1, and the integral in Eq. [3] thus saturates at the high-z limit. This is 
starkly different from the more intuitive continuous increase of optical depth with z, e.g., 
for Compton scattering and line absorption. The likely decrease of metallicity with z makes 
t igm(E,z) saturate even faster with z. The behavior of r IGM (E,z) without metallicity 
evolution and with r](z) oc (1 + z)~ 2 is shown in Fig. HJ overlaid on the measured opacities. 
These curves are not fits to the data, but are plotted to give a rough idea of how the IGM 
contribution to the opacity saturates at high z. 

The mean density of hydrogen in the IGM is n w 0.67 0^, (3 H 2 /87rGmfr) = 1.7 xl0~ 7 
cm -3 , where fi& = 0.045 is the Universe's baryon energy fraction (IKomatsu et al.ll2010l ). G is 



the gravitational constant, and m p the hydrogen mass. The pre-factor of 0.67 comes from the 
fact that ~ 74% of the Unive rse mass is in hydrogen atoms ~90% of which reside in the IGM 



( jFukugita and Peebles! 120041 ) . The Hubble constant is H = 71kms 1 Mpc 1 and the cross 



section at 0.5 keV <r(0.5 keV,0) = 6.22 x 10 22 cm 2 . For i](z) = 1 (no metall icity evolution). 



and u sing the standard cosmological parameters Qm = 0.27 and = 0.73 (jKomatsu et al. 



20101 ) . the integral in Eq. [3] up to z = 10 attains a value of ~1.4. It is interesting that this 
estimate results in the simple asymptotic expression tigm (0.5 keV, z ^> 1) ~ 2Z . In other 
words, a reasonable z = metallicity of Z ~ 0.2 could explain the measured asymptotic 
value of t ~ 0.4. 

There are two strong, oversimplifying assumptions in this estimate that have to do with 
the metallicity and the ion ization of the IGM . First, there is most definitely a metallicity 



evolution in the IGM (e.g., ISimcoe et al.l 120041 ) . If one assumes a universal metallicity evo 



lution of Z Q (z) = Z (l + z)~ k , the saturation of r IGM (E,z) occurs even at lower z than 
without evolution. In that case, for z' 3> 1, dr/dz' — > (1 + z')~ 2 ~ k , and the integral in Eq. [3] 
is approximately (1+k) times smaller. In other words, the corresponding expression is 
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r 7GM (0.5 keV, z > 1) w 2Z /(1 + fc) 



(4) 



For example, if one adopts for the IGM the metallicity trend of Z « (0.54 ± 0.10)(1 + 
Zq pa 0.5 and fc ~ 1.-25) ob served for Fe in the X-ray emitting gas of galaxy 



Z)-1.25± 0.25 (i e 

clusters up to z ~ 1 ( IBalestra et al.l 120071 ). the above relation suggests T/gm(0.5 keV, 2; ^> 
1) 0.45 ± 0.15, which is consistent with both the GRB mean value (above) and with 
the quasar results (next section). Alternatively, Z Q « (0.2 ± 0.1) (1 + z)' 1 derived for Fe 
from the measured supernova rate (Graur et al. 2011, private communications) suggests 
7"K?Af(0.5keV, z 3> 1) ~ 0.2 ±0.1, slightly less than the clusters or the above estimated GRB 
values. 

The second caveat is the ionization correction, which is more difficult to deal with, as the 
ioniz ation state of the diffuse IG M is not directly observed and is a matter of ongoing debate 
(e.g., iBolton Haehneltl 120071 ) . The above estimates are all based on the photo- ionization 



cross section of a neutral absorber. An ionized absorber would have a somewhat lower cross 
section at X-ray energies, and therefore a larger column density would be required to produce 
the same r. Calculating cross sections for different ionization states and different metallicities 
in the IGM is beyond the scope of this paper. However, in Fig. [1] we plot separately the 
contributions to the photo-ionization cross section of H and He, and that of the metals with 
solar composition. It can be seen that the metals are dominated by H and He below the O 
edge at 0.54 keV where they overtake the H and He contribution. For redshifted absorption 
it is thus the metal contribution that mostly determines r(0.5 keV) unless the metallicity is 
radically sub-solar. The pure metal contribution can also be thought of as an upper limit 
to the cross section for ionized plasma in which H and He are totally ionized (although in 
reality the metal contribution also somewhat decreases with ionization). A z = 1 absorber 
in which H and He are totally transparent (ionized) would still retain more than 70% of its 
opacity at 0.5 keV (observed) mostly due to C and O, and that fraction of course needs to 
be further scaled with the metallicity. 

To summarize this section, the mean high- z opacity for the diffuse IGM within the 
standard cosmological model can explain the soft X-ray opacity measured from spectra of 
high z GRBs (and quasars) as long as the gas is not too highly ionized, and as long as the 
metallicity is a reasonable fraction of the solar value and does not decrease too quickly with 
redshift. The low ionization is possible since the truly diffuse IGM suffered less gravitational 
collapse than the denser line absorbing IGM systems that have been heated to ~ 10 6 K 
( IDave et al.l 120011 ). The metallicity of 0.2 - 0.4 solar required to explain the absorption 
is actua lly lower than the Fe abundance observed in the hot gas of galaxy clusters up to 
z = 1.3 ( IBalestra et al.ll2007t ; iMaughan et al.ll2008l ). but slightly higher than that implied by 
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supernova rates (Graur et al. 2011) and th e abundances observed in IGM filaments towards 
nearby z < 0.4 metal absorbers (Z Q ~ 0.1, iDanforth fc Shullll2008l ). Recall that most of the 
IGM column is accumulated up to z = 1—2 (Eq.|2J Fig.H]), so even a sharp drop in metallicity 
beyond those redshifts does not change our con clusions. In fact, the mean metal licity in 
damped Lyman-a (DLA) absorbers up to z < 4 dSavaglioll2009l ; iKaplan etlUboioh is also 



consistent with these values, although the scatter in metallicit y in PL As (Prochaska et al. 



20031 ; ISavaglioll2009l ; iKaplan et al.ll2010l ). and in galaxy clusters (IBalestra et al.ll2007l ). at any 



given z is quite large. 



3.1. Host absorption 



It is impossible to rule out a cosmological evolutio n of the host (GRB en vironment or 
galaxy) column density that follows the values found by ICampana et al.l ( 120101 ) and reaches 



Nh ~ 10 cm" . On the other hand, at high redshift it is reasonable to assume that the 
contribution of the host galaxy to the opacity at 0.5 keV is small because of the redshift effect 
and the photo- ionization cross-section energy dependence a(E) ~ E~ 2 ' 5 discussed above. At 
low redshift the opacity of th e host certainly i s dom inant, but suffers from a large spread 
as is seen in Fig. |2l see also ICampana et al.l ( 120101 ) . This spread impedes any concrete 
conclusions regarding the hosts based only on their X-ray absorption properties. Such a 
spread is expected from the variety of host galaxies, GRB environments, GRB locations 
within the galaxies, and lines of sights. In Fig. H] we include a (high) mean host absorber 
of Nh = 3 x 10 21 cm -2 , whose contribution to r(0.5 keV) is significant at low redshift, but 
becomes less important than the IGM at z ~ 1.5. At high z, where the host contribution 
gradually becomes negligible, the spread in opacity values indeed becomes smaller, and the 
theoretical IGM opacity curve, thus, well describes the observed values. Since even for high-z 
sourc es, much of the IGM opacity is due to the gas up to z < 2 (see Eq. [3]), a clumpy IGM 
(e.g., 



Wiersma et al.l 1201 ll ) can also cause a spread in opacity at high z. 



4. Observed-frame opacity of X-ray selected high- z quasars 

We wish to further explore whether soft X-ray absorption of high-;? GRBs is due to 
the peculiar environment of the GRBs, or whether there is a significant intervening IGM 
contribution. A natural place to check for IGM absorption is with high-z steady X-ray 
sources, namely quasars. For this, we searched the XMM-Newton archive for all quasars with 
z > 2 for which the EPIC PN camera recorded more than 5000 source photons, to ensure 
high S/N. We found the 13 quasars listed in Tabled! most of which were already individually 
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reported in the literature. We used the standard EPIC PN pipeline products and fitted all 
the quasar spectra in the exact same manner as the GRB spectra, namely with a host (z) 
absorber. The column density Nh(z) obtained from this measurement was then translated 
into the optical depth r(0.5 keV) accor ding to Eq. [ TJ No a bundance correction was needed 
(Eq. |5J) as we used the abundances of Wilms et al.l ( 2000 ) for both the Galactic and host 
absorbers. The best-fit r(0.5 keV) values are listed in Tabled] and plotted in Fig. [Hand can 
be seen to compare well with the mean GRB opacities. Quoted and plotted errors correspond 
to the 90% confidence limits. When more than one XMM-Newton observation is available, 
we fit all observations simultaneously with a single absorbing column. 

The simple model of an absorbed (Galactic and host) power law fits most of the sources 
well as can be seen by the reduced x 2 /dof values quoted in Table [U and which are mostly 
close to unity. Three quasars in our sample, however, did not yield a satisfactory fit, namely 
QSOB0014+810, QSOB0438-43, and PKS 1830-210. The first two suffer fr om low cou nt 
statistics at high energies, but are fitted well up to E < 3 keV. C-statistic fits (]Cashlll979l ) of 
the entire spectra of these two sources are good and yield spectral parameters that are well 
within the errors of the limited band \ 2 fitting. We could not satisfactorily fit th e complex 
spectrum of the gravit at tonally lensed quasar PKS 1830-210 (IZhang et al.l 120081 ) with the 
same simple model, not even when the high energy region was ignored. Henceforth, we 
refer to our sample as the 12 targets in Table (H excluding PKS 1830-210. The sense of the 
absorption effect in this sample can be obtained from Fig. [5] In this figure, we plot the 
data to best-fit-model ratios for all targets of Table [T] except PKS 1830-210, but after the 
excess column has been removed from the model, i.e. Nh(z) is set to zero. When a source 
was observed more than once, we plot each spectrum separately (in color in the electronic 
version). The absorption effect can be seen to be consistent between different observations. 

For the most part, the optical depth toward quasars is well within the scatter of the GRB 
optical depths (Fig. [2]) and consistent with its mean trend (Fig. H]). However, the four lowest 
redshift targets with 2 < z < 2.5 do not show significant absorption. Three out of these four, 
essentially, have upper limits of r(0.5keV) < 0.1 (Table [T]). Formally, the quasar absorption 
effect is relatively uniform with an error-weighted mean of (r(0.5keV)) z> 2 = 0.19 ± 0.01, 
0.35 ± 0.09 with no weights, and 0.32 ± 0.08 if the one upper limit is also included, all of 
which are lower than the respective GRB values (Sec. [2]). The error-weighted mean optical 
depth of the eight quasars with z > 2.5 is (r(0.5keV)) 2> 2.5 = 0.41 ±0.02, which is consistent 
with the GRB value, and 0.43 ± 0.10 with no wei ghts. RBS315 (z = 2.69) has somewhat 



high absorption (see also iPiconcelli fc Guainazzil 120051 ) with respect to the quasars with 
r(0.5keV) = 1.08 ± 0.06, although this too is still well within the GRB scatter. It has been 
suggested that per haps the spectral curv ature in RBS315 is intrinsic to the source and not 
due to absorption (jTavecchio et al.ll2007l ). However, no definitive conclusion could be drawn. 
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The emergence of absorption at z > 2.5 raises an interesting connection with the re- 
i onization of He I I in the IGM observed from UV spectroscopy to occur at z = 2.7 ± 0.2 
( IShull et al.ll2010l ). The contribution of He to the observed opacity at 0.3 - 1 keV (1-4 keV 
in the rest frame for z = 2 — 3 targets) is small (Fig. [TJ, while the contribution of heavier 
elements such as C and O is appreciable. The ionization energy of He II is 54.4 eV, which 
is comparable to that of L-shell C and O charge states. For example, the ionization energy 
of O III is 54.9 eV, and that of C III is 47.9 eV. Thus, these elements are expected to be 
ionized into their L-shell when He II is ionized. The L-shell ionization of C and O results in 
their K-edge (responsible for the X-ray absorption) being pushed to higher energies, and in a 
slightly reduced photo-ion ization cross s e ction in the sub-keV regime. This was demonstrated 
recently for N in Fig. 8 of iGarcia et al.l ( 120091 ). The IGM X-ray opacity due to these metals 
will decrease more drastically once they become ionized down to their H-like state. This 
occurs for C and O at much higher ionization energies of 392.1 eV and 739.3 eV, respectively. 
If the reduced absorption observed along the lines of sight to the quasars at 2 < z < 2.5 
is due to the intervening material being highly ionized, it would conflict with ubiquitous 
IGM absorption, which is expected to nearly saturate by z ps 2.5 (see Sec. |3] and Fig. @|. 
Admittedly, the current quasar sample is small. A larger sample is needed to conclusively 
determine whether indeed X-ray absorption of high- z quasars increases significantly around 
z = 2.5. 



4.1. Radio quiet quasars 



Although the above sample was selected based on the availability of high S/N X-ray 
spectra in the XMM-Newton archive, the radio brightness of these bright X-ray sources creates 
a strong bias in the sample towards extremely radio loud quasars (RLQs), which are only 
a small fraction of the total quasar population. Are RQQs absorbed in the soft X-rays as 
much as the RLQs and GRBs? The answer to this question can not benefit from the high 
S/N available for the RLQs. Nevertheless, it is of considerable importance fo r assessing IGM 



absorption. Comp rehensive studies of high-z RQQs have been carried out by lShemmer et al. 



( 120051 . 120061 . 120081 ) . who find the observed soft X-ray spectrum of RQQs in these studies to be 
essentially featureless. In other words, significant absorption (or soft excess) is not detected. 
We obtained from the XMM-Newton archive the highest S/N spectra for RQQs. The data 
reductio n and handling i s sim i lar to that of the RLQs in Sec. HI We find similar results to 
those of IShemmer et al.l ( 120051 . 12006k 120081 ) and include these four data points in Fig. SJ 



In Fig. [6] we show the 4 RQQs with the most EPIC PN counts. The parameters of these 
sources are given in Table [2] and plotted in Fig. HJ It can be seen that the RQQ spectra 
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are much noisier than those of the RLQs and that none of the sources reveals significant 
absorption. The two RQQs with the most counts (bottom of Fig. [6]) are also at z < 2.5, 
where the RLQs do not show any significant absorption either (Figs. IU [5]). The two other 
RQQs at z > 4 suffer from even poorer S/N. Therefore, their spectra constrain the optical 
depth at 0.5 keV only poorly. The mean optical depth of ~ 0.3 found for the RLQs is still 
within the errors for these two sources, although their best fit values are lower. 

Despite the limited S/N compared to the R LQs (and many GRBs), the lack of absorption 
in RQQs (11 sources in IShemmer et al.ll2008l ) appears to be in contrast with the (diffuse 
IGM?) absorption effect observed in GRBs and RLQs. RQQs in the local universe are 
known to have soft X-ray excess emission that can offset the absorption effect. However, this 
excess that reaches ~ 1 keV in the quasar rest frame is not expected to affect the observed 
0.5 keV opacity for sources with z > 2, unless its strength increases with z. To summarize 
this section, the sense of no abs orption in RQQs is visible in Fig. [6] (see also the plots in 



Shemmer et al.ll2005l 12006k 120081 ). although some of the opacity upper limits are still high. 



5. Discussion and Conclusions 

The extragalactic soft X-ray opacity towards GRBs and quasars at low z is dominated 
by absorption in the host galaxy. However, the contribution of a fixed host column to the 
optical depth r at a given observed energy sharply drops with z, and the extragalactic opacity 
observed in the soft X-ray spectra of high-;? GRBs may thus have a dominant contribution 
from the diffuse IGM. Due to the redshift and energy dependence of the photo-electric cross 
section, the diffuse IGM opacity saturates at z > 2 and is expected to be isotropic and 
relatively constant beyond that redshift. The redshift where saturation occurs depends on 
the IGM metallicity. The IGM opacity is not expected to directly correlate with the line 
absorption systems observed in the optical and UV that originate either in the host galaxy 
or in over-dense IGM clumps. A large sample of Swift GRBs and a high S/N sample of RLQs 
observed with XMM-Newton appear to mostly be consistent with this picture of diffuse IGM 
absorption. However, a few RLQs at z < 2.5 and a few low S/N RQQs with very little 
absorption that nonetheless could represent a much larger quasar sample raise doubts. 

For the IGM to produce the X-ray opacities observed towards high- 2 GRBs and RLQs, 
it needs to be relatively enriched with metals and not too highly ionized. If indeed it is 
the diffuse IGM responsible for this absorption, it can account for a significant fracti on of 



the currently missing baryons, implied by big bang nucleosynthesis (jSteigmanl 120071 ). the 



observed angular power spectrum of the cosmi c microwave backgrou nd radiation and the 



Thomson opacity inferred from its polarization (IKomatsu et al.ll2010l ). Of these baryons, in 
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the local universe only ~ 50% are present in the galaxies, galaxy clusters and UV-optical 
IGM line systems known to date (for a review see iBregmanl 120071 ) . Furthermore, if 90% of 
the baryons are in the IGM and if the diffuse low-z IGM metallicity is indeed ~ 0.2 — 0.4 
solar as we postulate here, then the IGM contains the bulk of the metals in the present-day 
universe, compared to, say, solar metallicity in stars and galaxies that comprise only 10% of 
the baryons. 

The presence of absorption in GRBs and RLQs, that both harbor powerful jets, but 
less in RQQs, also raises the possibili ty that the absorption effect has so mething to do with 
the jet. It was discovered by Rosat fab/is et al. 1994 ; Fiore et al. 1998 ) that high-z RLQs 
are much more absorbed in the X-rays than RQQs, which seems to support a jet effect. On 
the other hand, these authors also realized the absorption (if intrinsic) increases with z and 
not with luminosity, which argues against a jet-physics origin. The fact that local RLQs 
generally do not have the high column d e nsitie s that are found in the high-z sources has 
been confirmed recently by iGalbiati et al.l (120051 . Table 4 therein). Indeed, ascribing photo- 
electric absorption to the jet is counter-intuitive, as the jet is not expected to comprise atomic 
material with bound electrons, and especially not metals, that produce the observed X-ray 
opacity. Moreover, the little scatter of the observed optical depth found in this paper, and 
even more the tendency to an asymptotic opacity at high z that require a putative intrinsic 
column to scale approximately with (1 + z) 2 ' 5 to offset the decreasing cross section, calls into 
question the realistic role jets can play in determining the soft X-ray opacity. 

Turning to the trends of intergalactic line absorption towards GRBs and quasars in the 
optical band does not provide a clear-cut answer to the uniqueness of absorption towards 
jetted sources either. On one hand , the numbe r of MgH intervenin g absorption systems 
towards GRBs ( Prochter et al.ll2006l ) and blazars ([Bergeron et abj|2010j ) is a few times higher 
than that towards RQQs. As there is no obvious reason for line sights towards GRBs to 
be different from those towards quasars, it would require metals to be entrained in the jets 
of GRBs and RLQs. On the other hand, line sights towards high-z RLQs seem to have 
a similar number density of in tervening DLA syst ems as those towards optically-selected 
(again, mostly RQQs) samples (jEUision et al.l 1200 ll ). 



Prospectively, there is a clear way to distinguish a well confined absorber from a cos- 
mologically diffuse one based on line absorption. The detection of absorption lines, or lack 
thereof at high spectral resolution (and high S/N) would provide a definitive characterization 
of the absorber. The CCD X-ray detectors used in this work do not have the spectral resolu- 
tion required to discern lines, while the grating spectrometers on board XMM-Newton do not 
have the sufficient effective area to provide adequate spectra, at least not with the modest 
exposures available to date. High S/N high- resolution spectra of the quasars either through 
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a long XMM-Newton exposure, or with future instruments should be able to unambiguously 
determine whether X-ray absorption of high- ,2 sources is due to their hosts or due to the 
diffuse IGM. 

The authors thank A. Heller for assistance with the GRB data extraction, and S. Kaspi 
and A. Nusser for early discussions. We acknowledge important insights by D. Chelouche on 
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Table 1. High- z Quasars 



Source 


XMM-Newton 


z 


photons 


N H (G&\.) 


N„(z) a 


t(0.5 keV) b 


X 2 /dof 




Obs ID(s) 






10 20 cm- 2 


10 22 cm- 2 










CjJdd rS14zc5+4z± i 


L)111zdU1U1/7U1 


4.715 


12580 


1 1 O 


n i i n rr 
Z.l ± U.O 


0.24 


± 


0.06 


1 n a 
1.U4 




UZlZ4oU iUl 


















T3A/TAT TO^OK QQ/IQ 
r^lVllN JU0z0-oo4o 


UUoUloUlUl/ oUl 


A A 1 Q 


ooonn 
zooUU 


L.LO 


1 o j- n q 
i.y ± \j.o 


0.26 


± 


0.04 


U.yo 




m /i nrcnn^m /7ni /om /nm 

ui4yoUL)DUi/ /ui/oUi/yui 




















U14yoUlUUl/ 1U1 


















RXJ 1028.6-0844 


0093160701, 0153290101 


4.276 


6250 


4.49 


1.8 ± 0.7 


0.27 


± 


0.10 


0.98 


QSOB0014+810 c 


0112620201 


3.366 


12500 


13.6 


1.8 ± 0.5 


0.41 


dz 


0.11 


0.98 


PKS 2126-158 


0103060101 


3.268 


35200 


4.92 


1.6 ± 0.2 


0.39 


± 


0.05 


0.98 


QSOB0438-43 c 


0104860201 


2.852 


7400 


1.35 


1.8 ± 0.3 


0.59 


dz 


0.12 


1.08 


RBS315 


0150180101 


2.690 


69950 


9.26 


2.9 ± 0.2 


1.08 


± 


0.06 


1.00 


PKS 2351-154 


0203240201 


2.675 


9000 


2.51 


0.6 ± 0.2 


0.21 


dz 


0.07 


1.05 


PKS 1830-210 


0204580201/301/401 


2.507 


75500 


20.2 


18 ± 0.6 


7.1 


dz 


0.2 


1.18 


QSO J0555+3948 


0300630101 


2.363 


5000 


28.2 


0.49 ± 0.44 


0.22 


± 


0.20 


0.97 


PKS 2149-306 


0103060401 


2.345 


36200 


1.61 


0.08 ± 0.07 


0.036 


dz 


0.032 


1.00 


PKS 0237-230 


0300630301 


2.225 


12550 


2.16 


0.10 ± 0.10 


0.05 


± 


0.05 


1.06 


4C 71.07 


0112620101 


2.172 


225450 


2.85 


0.09 ± 0.03 


0.05 


dz 


0.02 


1.08 



a Measured assuming an absorber at the host at redshift z. 
b Deduced from Nh(z) using Eq. [TJ 
c Fitted only up to 3 keV. 



Table 2. Radio Quiet Quasars 



Source XMM-Newton z photons N H (G&\.) N H (z) a r(0.5 keV) b x 2 /dof 





Obs ID(s) 






10 20 cm- 2 


10 22 cm- 2 






PSS 0926+3055 
Q 0000-263 
HE 2217-2818 
Q 1318-113 


0200730101 
0103060301 
0302380401 
0402070301 


4.190 
4.111 
2.414 
2.306 


800 
1100 
2250 
2500 


1.89 
1.67 
1.28 
2.22 


n ,: !5 +1 - 5 

u -°°-0.35 
5+ 13 

u -°-0.5 

+o.i 

fi 1+0.26 


0.06±S5S 

n +0.05 
u -0 

05+ 012 


1.02 
1.26 
0.90 
1.07 



a Measured assuming an absorber at the host at redshift z 
b Deduced from Nh{z) using Eq. [T] 
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0.1 1 10 



Photon Energy (keV) 



Fig. 1. — Photo-ioniza t ion cr oss section per H atom extracted from the Xspec phabs model 
using the IWilms et all (120001 ) abundances preferred in this work (solid line), and separated 
into its H and He contri bution (dashed) and metal contribution (dash-dot). The cross 
section obtained with the lAnders and Grevessd (119891 ) abundances (dotted line) are shown 
for comparison. 
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Fig. 2. — Optical depth r at 0.5 keV as a function of redshift for the Swift/XRT PC-mean 
GRB sample. Black squares are detections and red stars represent error-weighted (r) values 
averaged over Az = 1 bins; Note the striking constancy of (r) with z and its tendency 
towards the value of ~ 0.4 for z > 2. Blue triangles represent +90% confidence (upper) 
limits. 



- 19 - 



s i - 

° «! rr 



GRB 090423 
z=8.26 972 cts 



I I 



ft 



GRB 050904 
z=6.29 5096 cts 



- + 



• it 



iU.. 

lit" 1 ' 



GRB 060522 
z=5.11 564 cts 



3* 



GRB 06051 OB 
■ z=4.9 462 cts 






i 







GRB 100513A 
z=4.77 523 cts 



GRB 060223A 
z=4.41 391 cts 



I- 



+ 





0.5 1 2 5 

Observed Energy (keV) 



1 2 5 

Observed Energy (keV) 



1 2 5 

Observed Energy (keV) 



Fig. 3. — Data to model ratio plots for the twelve highest-z GRBs with confirmed absorption. 
Data are binned to conveniently represent the extra-galactic transmission functions. Note 
the overall similar absorption amplitude irrespective of z. 
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1+Z 

Fig. 4. — Optical depth r at 0.5 keV as a function of redshift. Mean Swift/XKT GRB values 
are denoted by red stars (c.f., Fig. [2]). Black (green) circles (triangles) represent the RLQs 
(RQQs) measured with XMM-Newton. Upper limits are denoted by down pointing triangles 
(RQQs are all upper limits). Theoretical IGM and host-galaxy contributions (Sec. [3]) to 
the total opacity (black line) are shown separately. The IGM theoretical contribution is 
scaled to approach r = 0.4 at high z. The dashed line represents the IGM contribution if its 
metallicity evolved as Z & (z) = Zq(1 + z)~ 2 . 
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Fig. 5. — Data to model ratio plots for the QSO sample of Table [TJ Data are binned to 
conveniently represent the extra-galactic transmission functions. Note the overall similar 
absorption effect, but the lack of absorption for z < 2.5. Multiple spectra for a given source 
represent separate XMM-Newton observations. 



- 22 - 



- 1 1 1 1 — i — r- 



PSS 0926+3055 
z=4.19 792 cts 



CO 

CC 
"55 

"O ,_ 

o 



B 

Q in 
o 



+" 



- 1 1 — i — i — i — r- 



- 1 1 1 — i — i — i — r- 



Q0000-263 
z=4.10 1079 cts 



0.5 1 

— I 1 1 1 1 — I — | 

Q1318-113 
z=2.31 2504 cts 



5 1 

— I 1 1 1 — I — 




4 



++ 



0.5 1 2 

Observed Energy (keV) 



i i i i i i i 

0.5 1 



—I i I i i i l_ 



10 



-I 1 1 1 1 — I — | 

HE2217-2818 
z=2.41 2229 cts 



-i 1 1 — i — i — r 



: 'l' "Wj'l fjf't 



10 



0.5 1 2 

Observed Energy (keV) 



10 



Fig. 6. — Data to model ratio plots for the highest-count RQQs listed in Table [2] Plotted are 
in effect the extra-galactic transmission functions. Data are binned to conveniently represent 
the extra-galactic transmission functions. The lack of absorption seen for the two z < 2.5 
sources is consistent with the RLQ results (Fig. [5]), while for z > 4 it is not. This result, 
however, needs to be taken with caution due to the low S/N RQQ spectra. 



